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1. INTRODUCTION 

Cardiovascular disease is a collection of irregularities affecting the heart; it is considered one of the 
most important causes of death in the world. According to the World Health Organization, there were 
approximately 17.9 million deaths in 2016; this big number represents 31% of deaths worldwide [1]. The 
lives of people with cardiovascular diseases are in constant danger, quick and effective diagnosis of these 
diseases can save a lot of lives. Several techniques in the medical field are used to diagnose cardiovascular 
disease, such as blood tests, coronary angiography, cardiac MRI, X-ray and electrocardiography. However, 
most of these techniques require medical assistance from experienced people, which is not always the case if 
we knew that almost 30% of cases with these diseases come from poor countries. Electrocardiography is a 
non-invasive detection technique based on recording the electrical activity of the heart over time [2], the 
signal obtained during this recording is called an electrocardiogram (ECG). The ECG signal is considered 
among the most widely used biomedical signals to detect heart problems, ECG signal contains a large 
number of information that can be of great interest in the detection and diagnosis of many heart diseases [3] 
which appears in some distortions of the signal shape. Despite the advantages of the use of ECG signals, 
there are many limitations of this technique; the difficulty of interpreting the signal and the lack of 
experienced personnel are among the constraints most encountered during identification of ECG signals; 
also, ECG signal contains various unwanted noises that prevent the correct extraction of useful and necessary 
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information for classification [2], [4]. Hence, finding solutions to solve these problems has become a 
necessity. The field of medical engineering try the challenge by developing models capable of processing 
ECG signals and dentifying any abnormalities present in the signal [5], thus, signal processing is required; 
signal processing techniques analyze efficiently various kinds of signals, especially ECG signals. The aim of ECG 
signal processing is to extract features to distinguish between normal signals and those representing abnormalities. 

Several techniques based on signal processing have been developed over years for processing ECG 
signals. Ahlstrom and Tompkins have used digital filters for real time ECG signal processing [6] to denoise 
the signal and detect the QRS complex. Hargittai [7] used Savitzky-Golay Least-Squares Polynomial filters 
to preserve the details of the signal, Francisco et al. [8] process ECG signal using principal component 
analysis (PCA), Haque et al. [9] use adaptive filtering algorithms, while Ahmed et al. [10] use a method 
based on cross correlation theory. Gustavo et al. [11] propose a comparison between several methods to 
remove baseline wander. Also the choice of features and the method used for its extraction affect directly the 
quality of signal characterization, some techniques are based on the extraction of morphological features such 
as the detection of the QRS complex proposed by Jiapu [12], the calculation of R-R intervals and peaks 
detection mentioned in the works of Shanti Chandra et al. [13] and Priyanka [14], other methods are based on 
the extraction of the statistical ones or even a mixture of DWT and statistical features, Abdullah et al. [15], or 
morphological and statistical features, Sahoo et al. [16]. 

The discrete wavelet transform (DWT) is a very powerful tool in the field of signal processing [17]- 
[20], this technique gives satisfying results in the processing of the noise which affects the signal that can 
allows us to reconstruct a denoised signal, also the DWT allows to extract different features that characterize 
the signal. The interest of the discrete wavelet transforms (DWT) pushes us to better exploit this technique in 
the processing of ECG signal [21], [22]. Also, machine learning techniques are constantly evolving; this 
evolution is reflected in the use of these techniques in several areas [23], more precisely in the classification 
and identification of signals. The combination of signal processing techniques and machine learning models 
gives us promising results and the performance of each model changes depending on the chosen algorithm 

The purpose of this modest work is to establish a characterization model of ECG signals able to 
differentiate between normal and abnormal signals, the adopted model is based on the extraction of statistical 
features from the approximation coefficients obtained by the wavelet decomposition of the signal and classify 
these features using an SVM classifier; this model was approved after the choice of the best wavelet and the 
best scale increasing it accuracy, the recordings of the used ECG signals are taken from the MIT-BIH 
arrhythmia database [24], [25]. It’s composed of 48 recordings of different patients; these patients are 
classified into two categories: A healthy and a sick one. After the processing of the ECG signals and the 
extraction of features a classification of patients to healthy or sick ones is needed. To do this we have chosen 
a support vector machine (SVM) classifier who seems most appropriate to do this task [26]-[28]. 


2. PRELIMINARIES 
2.1. Wavelet analysis 
Wavelet analysis is one of the powerful tools in signal processing, it is considered as a technique 
which aims to solve the problem of non-stationary signals. This notion was introduced in the 20th century by 
Haar who constructed the simplest wavelet, and then developed in the 1980s through the research work of 
Mallat [18], [19], [29]. One of the powerful things of this approach is that it allows a signal to be analyzed in 
time and frequency, which makes it very useful in extracting the various information contained in that signal. 
The wavelet transform decomposes a signal x(t) using a series of wavelets YW which derive from a 
mother wavelet by a displacement in time carried out by a translation and a dilation or compression by 
scaling of the mother wavelet, this series of wavelets is defined as [30]: 


apt) = =H) (1) 


where a and b are respectively the scaling and the translation coefficients W being the mother wavelet which 
must verify the following condition [17]: 


se pat = 0 (2) 


The wavelet transform exists in two principal forms the CWT and the discrete wavelet transform (DWT), 
thus we define the continuous wavelet transform (CWT) of a signal x(t) as [18]: 


Cop = ln x(OW" (Dat 3) 
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W* represents the complex conjugate of the wavelet W (t). 

The computation of this continuous wavelet transform seems sometimes difficult and time 
consuming. The discrete wavelet transform give a solution to this problem. To do this we use an algorithm 
(proposed by Mallat) [18] to calculate the discrete wavelet coefficients, this decomposition consists in 
passing the signal through a succession of complementary low-pass h[n] and high-pass g[n] filters (in stages), 
the low-frequency components resulting from the low-pass filters h[n] represent approximation coefficients 
and those from the high frequency g [n] filter are called the detail coefficients [31] as shown in Figure 1. 


2.2. Features extraction 
The choice of the nature, number and method of features extraction is considered one of the decisive 
steps in the characterization of a signal, since each of these parameters contain various information which 
influence on the efficiency of classification. In the literature many are the features and the methods used to 
characterize a signal especially for the ECG signals, there are the morphological features such as the PR, RR and 
QRS intervals, which are determined basically by the detection of the QRS complex [12]-[14]. There are also 
statistical features such as maximum, minimum, mean, variance and standard deviation [15], [16], [32], [33]. 
In this related work we will focus on the following features: 
The mean value (4): 


p= Thx; (4) 


~ oN 


The root mean square (5): 


1 
rms = ove ne (5) 


The variance (6): 
1 
v= =YN 05 — u)? (6) 
Standard deviation (7): 


gs Fates — py? (7) 


Kurtosis (8): 


NN (xj-p)3 
g= ee (8) 
Wavelet entropy (9): 
H = ~~ 5M, xPlog(x?) (9) 


where N is the length of the signal x 


Figure 1. Signal decomposition using DWT at 3rd level of scale 
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2.3. Classification 

The chosen classification technique is support vector machine (SVM), SVM is a machine learning 
technique that allows binary classifications. The purpose of this method is to separate the features of the 
signal into different categories [27], [28]. The purpose of this method is to separate the features of the signal 
into different categories, it's based on the construction of a separation area between the different classes of 
the learning set{x;,y;}, with y; € {-1,1} and i = 1,2,....,n the lines that delimit this area are called a 
hyperplane which is defined by [28], [29]: 


y=wxt+b (10) 


W represents a normal vector on the hyperplane and b is the bias. 


Thus, the area where we have: y; = —1 implies that wx; + b < —-1 (11) 
Similarly, for y; = +1 we have wx; +b = 1 (12) 
In all cases we have y;(wx; + b) = 1 (13) 
Good separation aims at increasing the width of the margin between hyperplanes, 

this is equivalent to minimizing a knowing thaty; (wx; + b) -—1 = 0. (14) 


If we can't separate linearly these classes, then the optimization problem becomes [26]: 


IIwIl? 


min(—— + € di-1 i) 
avec y; (wx, +b)-& =1 


(15) 


3. METHODOLOGY 
In this work we want to establish a simple model that allows: 
— ECG signal processing. 
— Extraction of the features that characterize the signal, these features will be calculated by applying 
(4)-(9). 
— The classification of signals processed into normal (healthy patients) and other abnormal signals (the 
sick ones) using a linear SVM classifier. 
Processing techniques and extraction methods are based on discrete wavelet analysis; Figure 2 
shows the process of identifying ECG signals. 


3.1. Database 

The ECG signals are taken from the MIT BIH data base of physionet [24], [25], this database 
contains 48 recordings each recording is half hour length, the recordings are sampled at 360 Hz, in this study 
we analyzed duration of one minute for each signal. 


3.2. ECG signal processing 

The recordings of ECG are not protected from noises; these noises cause significant perturbations 
during classification and diagnosis. For ECG signals, there are different sources of noises, noises of technical 
origins such powerline interference (50 or 60 Hz noise from mains supply) or bad wiring and others of 
physical origins generated by the physical activity of the body such electromyographic (EMG) noise, 
baseline wander (low frequencies) which are considered among the most important [34]. The steps of signal 
denoising are shown in Figure 3. 

In this stage of processing, we tried to eliminate these noises using the discrete wavelet 
decomposition. After the decomposition of the signal into eight levels we notice that the noise (high 
frequencies) is located in the detail coefficients d; and d2. For the baseline wander we notice that the 
approximation coefficient Ag corresponds to the frequency interval of these fluctuations 0-0.5 Hz [20]; 
Figure 4 shows detail and approximation coefficients resulting from the wavelet decomposition of recording 
121 m, and Figure 5 shows part of the same recording before, as shown in Figure 5(a) and after, as shown in 
Figure 5(b) denoising, the mother wavelet used in the decomposition is Symlet [35]. 
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Figure 2. Signal denoising process Figure 3. Diagram of the proposed method 
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Figure 4. Wavelet decomposition of ECG signal (record No 121) d1, d2, d3, d4, d5, d6, d7 and a8 


3.3. Features extraction 

After denoising the ECG signal, we decompose the denoised signal using DWT by using different 
types of wavelets to analyze the signal, and then the approximation coefficients are calculated up to scale 8 as 
shown in Figure 6. After the extraction of the approximation coefficients, the features mentioned above (4)- 
(8) are calculated from the approximation coefficients a1, a2, a3, a4, a5, a6, 47, as, thus we create for each level 
a data set containing the features for all the records. 


3.4. Classification 

The resulting dataset is introduced in a SVM classifier and it’s divided into two subsets; the first will 
be served for training the model, while the second will be designed to test the performance of the model. The 
choice of the training and test subsets is crucial and it can causes many problems which influence the 
effectiveness of the model such as overfitting and underfitting [36]; to limit this kinds of problems we use k 
flod cross validation to select this subsets [37], [38], this method consists in dividing the data set into parts of 
number of &, one of these parts is chosen as a subset of test and the other k — 1 parts as a training subset, this 
process is repeated k times and each time a different part is taken to do the test. The precision of the model is 
obtained by averaging the precision of each iteration [39]-[41]. In this work we took k = 8, that is to say that 
7/8 of the data set is taken for training and //8 is considered as a test subset. The next step will be devoted to 
measuring the performance of the model for the different types of mother wavelets at different scales (up to 
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the 8th scale); the accuracy (Acc), the sensitivity (Sen) and the specificity (Spe) are the quantities that inform 
us about the performance of each model with [29]: 


TN+TP 


Acc = —————-_ x 100 (16) 
TP+TN+FP+FN 
Sen = —"-x 100 (17) 
TP+FN 
TN 
Spe = —*— x 100 (18) 


TP: True positive (represents normal signal who were correctly classified) 

TN: True negative (represents abnormal signal who were correctly classified) 
FP: False positive (represents normal signal who were incorrectly classified) 
FN _: False negative (represents abnormal signal who were incorrectly classified) 
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time(s) 
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Figure 5. Record no 121 before and after denoising; (a) original signal and baseline wander, (b) signal after 
removing noise and baseline wander 


Electrocardiogram signals classification using discrete wavelet transform and... (Youssef Toulni) 


966 O ISSN: 2252-8938 


al a5 
400 1000 
200 600 
0 0 
200 -500 
0 2000 4000 Samp 8000 10000 12000 0 100 200 300 ab 400 500 600 700 
600 600 
ee 
600 -600 
0 1000 2000 30% 4000 6000 6000 0 60 100 160 a7 200 260 300 350 
4000 500 
500 
0 
0 
600 600 
0 500 1000 150) 2000 2600 3000 0 20 a) 60 80 38 100 120 140 160 #8 180 
1000 500 
600 0 
0 600 
500 -1000 
0 200 400 600 800 1000 1200 1400 0 10 20 30 a 50 60 70 80 90 


Samples Samples 


Figure 6. Approximation coefficients al, a2, a3, a4, a5, a6, a7, a8 


4. RESULTS 

Table 1 as shown in appendix, summarizes the results obtained for different wavelet families at the 
first eight scales of the approximation coefficients. As we can see, the best accuracy is generally obtained at 
the fourth scale. The study carried out by Siti [42] which is based on the calculation of statistical features 
from the detail coefficients at the 4th, 5th and 6th and the KNN as a classification method reaches an 
accuracy of 71% using the sym7 wavelet; while in this work we had reached an accuracy of 81.67% using 
only the statistical features extracted from the approximation coefficient a4 for the same wavelet and an SVM 
classifier; ;also in the study of Siti [42] we reached an accuracy of 85% as the best accuracy by calculating 
the MFCC coefficients, while in this work we found 87.50% as better accuracy by adopting only the coif5 
wavelet as the mother wavelet as shown in Table 2. 


Table 2. Comparison results 


Method Classifier _ Accuracy 
Haar 68% 
Siti [42] DE DEER DO coe TENN N% 
MFCC 85% 
Proposed method Sym7 81,67% 
POE cues OY _ersog 


5. CONCLUSION 

In this article we focused on the classification of ECG signals, we were then able to establish a 
model that allows to process and classify these signals. ECG signals were processed using DWT discrete 
wavelet analysis, the discrete wavelet transform allows the calculation of approximation coefficients which 
are used to extract the features for different scales, the dataset which consists of these parameters is fed into 
an SVM classifier with cross validation in order to distinguish between a normal and an abnormal signal. 
This model was tested with different mother wavelets and at different scales to reach an accuracy of 87.50%, 
which was the best. To conclude, the choice of the mother wavelet, the scale of the decomposition as well as 
the size of the training and test sets have a considerable influence on the accuracy of the model. 
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Table 1. Accuracy, sensitivity and specificity of the different wavelets at the first eight scales (continue) 


Wavelet Classifier Approximation coefficients scale level 
performance 1 Dy 3 4 5 6 7 8 
Accuracy 72.08 71.67 70.83 79.58 70.42 68.33 60.42 56.84 
db2 Sensitivity 75.83 76.67 77.50 90.00 81.67 75.42 68.33 64.50 
Specificity 68.33 66.67 64.17 69.17 59.17 60.00 52.50 50.00 
Accuracy 71.25 72.92 72.92 77.92 70.42 65.42 58.75 53.51 
db3 Sensitivity 74.17 76.67 79.17 86.67 79.17 74.58 69.17 53.83 
Specificity 68.33 69.17 66.67 69.17 61.67 53.33 48.33 53.33 
Accuracy 72.08 71.25 T4.17 75.83 68.33 64.17 58.75 51.97 
db4 Sensitivity 75.00 75.83 80.00 86.67 78.33 72.50 64.17 52.54 
Specificity 69.17 66.67 68.33 65.00 58.33 54.17 53.33 51.67 
Accuracy 72.92 70.83 71.25 77.08 68.33 65.00 60.00 39.99 
db5 Sensitivity 79.17 76.67 79.17 89.17 80.00 74.17 69.17 55.53 
Specificity 66.67 65.00 63.33 65.00 56.67 55.00 50.83 29.17 
Accuracy 73.33 72.50 72.92 76.25 65.00 65.83 59.58 41.08 
db6 Sensitivity 80.00 78.33 80.00 87.50 74.17 75.42 73.33 68.19 
Specificity 66.67 66.67 65.83 65.00 55.83 53.33 45.83 18.33 
Accuracy 70.83 73.33 77.50 72.08 67.08 63.33 60.00 41.68 
db7 Sensitivity 78.33 79.17 85.00 81.67 75.83 75.00 74.17 69.26 
Specificity 63.33 67.50 70.00 62.50 58.33 49.17 45.83 18.33 
Accuracy 73.33 74.58 73.75 75.83 67.50 63.75 59.58 49.58 
db8 Sensitivity 79.17 81.67 80.83 86.67 73.33 75.42 67.50 68.59 
Specificity 67.50 67.50 66.67 65.00 61.67 48.33 51.67 30.00 
Accuracy 72.92 71.67 72.92 77.50 70.00 67.92 60.83 54.49 
sym2 Sensitivity 76.67 75.00 77.50 85.00 80.00 76.25 71.67 65.10 
Specificity 69.17 68.33 68.33 70.00 60.00 57.50 50.00 45.83 
Accuracy 71.67 71.25 72.08 77.08 70.83 64.17 60.83 56.30 
sym3 Sensitivity 76.67 75.00 77.50 89.17 80.00 7292 72.50 58.28 
Specificity 66.67 67.50 66.67 65.00 61.67 52.50 49.17 55.00 
Accuracy 72.92 71.25 71.67 81.67 67.50 72.92 65.83 45.58 
sym4 Sensitivity 75.83 74.17 77.50 95.00 80.83 79.58 70.83 64.00 
Specificity 70.00 68.33 65.83 68.33 54.17 64.17 60.83 30.83 
Accuracy 70.00 71.67 72.50 82.08 71.25 72.08 63.33 44.36 
sym5 Sensitivity 75.83 75.00 75.83 93.33 85.83 81.25 68.33 59.91 
Specificity 64.17 68.33 69.17 70.83 56.67 60.83 58.33 31.67 
Accuracy 71.67 T4AAT 71.25 79.17 71.67 72.50 61.67 46.51 
sym6 Sensitivity 76.67 77.50 76.67 91.67 85.00 81.25 69.17 61.18 
Specificity 66.67 70.83 65.83 66.67 58.33 61.67 54.17 35.00 
Accuracy 71.67 72.50 72.50 81.67 72.08 70.83 58.75 42.00 
sym7 Sensitivity 76.67 76.67 75.83 91.67 85.83 81.25 66.67 58.77 
Specificity 66.67 68.33 69.17 71.67 58.33 57.50 50.83 27.50 
Accuracy 73.33 72.92 73.75 79.17 68.75 67.50 60.00 42.98 
sym8 Sensitivity 78.33 75.83 77.50 90.00 80.00 77.50 69.17 63.21 
Specificity 68.33 70.00 70.00 68.33 57.50 55.83 50.83 25.83 
Accuracy 72.08 72.50 71.67 79.17 71.67 72.92 66.67 51.72 
coifl Sensitivity 75.00 75.83 75.00 88.33 82.50 80.00 70.83 65.32 
Specificity 69.17 69.17 68.33 70.00 60.83 64.17 62.50 40.00 
Accuracy 72.08 72.50 71.67 79.17 71.67 72.92 66.67 51.72 
coif2 Sensitivity 75.00 75.83 75.00 88.33 82.50 80.00 70.83 65.32 
Specificity 69.17 69.17 68.33 70.00 60.83 64.17 62.50 40.00 
Accuracy 72.08 72.50 71.67 79.17 71.67 72.92 66.67 51.72 
coif3 Sensitivity 75.00 75.83 75.00 88.33 82.50 80.00 70.83 65.32 
Specificity 69.17 69.17 68.33 70.00 60.83 64.17 62.50 40.00 
Accuracy 74.48 73.96 75.00 85.42 71.88 72.40 60.42 49.40 
coif4 Sensitivity 81.25 81.25 83.33 97.92 83.33 79.69 72.92 63.96 
Specificity 67.71 66.67 66.67 72.92 60.42 62.50 47.92 38.54 
Accuracy 75.00 73.96 75.00 8750 73.96 73.96 55.21 42.80 
coifS Sensitivity 82.29 80.21 82.29 100.00 84.38 82.29 65.63 55.31 
Specificity 67.71 67.71 67.71 75.00 63.54 63.54 44.79 33.33 
Accuracy 63.75 65.83 67.50 70.42 69.17 64.58 62.92 61.70 
biorl.1 Sensitivity 68.33 69.17 72.50 79.17 85.00 72.92 75.83 73.63 
Specificity 59.17 62.50 62.50 61.67 53.33 53.33 50.00 50.83 
Accuracy 71.67 71.25 73.33 79.17 69.58 70.42 60.83 47.44 
bior1.3 Sensitivity 75.00 75.00 76.67 89.17 83.33 76.25 70.00 62.57 
Specificity 68.33 67.50 70.00 69.17 55.83 63.33 51.67 35.00 
Accuracy 70.42 71.25 71.25 78.75 72.50 74.17 60.42 40.99 
bior1.5 Sensitivity 74.17 75.00 75.83 90.00 85.83 84.17 66.67 61.32 
Specificity 66.67 67.50 66.67 67.50 59.17 61.67 54.17 22.50 
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Table 1. Accuracy, sensitivity and specificity of the different wavelets at the first eight scales 
Classifier Approximation coefficients scale level 
performance 1 1 1 
Accuracy 72.08 72.08 69.58 72.50 77.08 69.58 66.67 54.28 
bior2.2 Sensitivity 75.00 75.00 T4AAT 80.00 89.17 75.83 64.17 61.56 
Specificity 69.17 69.17 65.00 65.00 65.00 61.67 69.17 47.50 
Accuracy 71.67 70.83 70.42 77.50 70.42 72.08 60.83 49.27 
bior2.4 Sensitivity 76.67 T4AAT T4AAT 85.83 80.83 76.25 69.17 62.49 
Specificity 66.67 67.50 66.67 69.17 60.00 66.67 52.50 38.33 
Accuracy 70.00 71.25 70.42 7917 72.50 T4AAT 61.25 46.13 
bior2.6 Sensitivity 75.00 75.83 75.83 88.33 83.33 80.83 67.50 57.49 
Specificity 65.00 66.67 65.00 70.00 61.67 65.83 55.00 37.50 
Accuracy 70.42 71.67 70.42 78.33 T3315 72.92 62.08 43.62 
bior2.8 Sensitivity 75.00 75.83 72.50 86.67 85.00 81.25 70.00 64.06 
Specificity 65.83 67.50 68.33 70.00 62.50 62.50 54.17 26.67 
Accuracy 63.33 64.17 65.83 69.17 70.42 62.50 64.58 60.58 
rbiol.1 Sensitivity 67.50 67.50 69.17 78.33 85.00 72.50 78.33 76.17 
Specificity 59.17 60.83 62.50 60.00 55.83 49.17 50.83 46.67 
Accuracy 70.00 72.50 72.92 77.08 72.92 70.42 66.67 60.36 
rbiol.3 Sensitivity 75.00 75.00 78.33 85.00 82.50 77.50 80.83 75.53 
Specificity 65.00 70.00 67.50 69.17 63.33 61.67 52.50 48.33 
Accuracy 69.17 71.67 72.92 78.33 T4AAT 72.08 64.17 52.65 
rbiol.5 Sensitivity 73.33 75.83 78.33 87.50 81.67 80.42 75.83 70.90 
Specificity 65.00 67.50 67.50 69.17 66.67 61.67 52.50 37.50 
Accuracy 72.92 72.50 75.83 76.67 70.42 70.00 57.08 52.33 
rbio2.2 Sensitivity 75.83 77.50 84.17 85.83 81.67 80.83 62.50 54.11 
Specificity 70.00 67.50 67.50 67.50 59.17 56.67 51.67 50.83 
Accuracy 71.67 72.08 75.83 78.33 68.33 68.75 56.67 44.21 
rbio2.4 Sensitivity 76.67 77.50 83.33 87.50 78.33 73.33 66.67 56.94 
Specificity 66.67 66.67 68.33 69.17 58.33 62.50 46.67 34.17 
Accuracy 70.42 71.25 77.50 75.42 70.00 68.75 54.17 40.63 
rbi02.6 Sensitivity 75.83 76.67 85.83 83.33 TIAT 75.83 65.83 56.03 
Specificity 65.00 65.83 69.17 67.50 60.83 60.83 42.50 28.33 
Accuracy 72.92 T4AAT 78.33 7917 72.92 65.83 50.83 45.07 
rbio2.8 Sensitivity 77.50 81.67 85.83 90.83 81.67 73.33 65.83 61.33 
Specificity 68.33 66.67 70.83 67.50 64.17 55.83 35.83 30.83 
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